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SUMMARY 


A  mathematical  model  of  the  Red  River  of  the  North  Basin, 
Minnesota  and  North  Dakota,  can  be  used  as  a  water  quality  planning 
management  tool  to  simulate  time  and  spatial  variations  of  flovr  and 
concentrations  of  total  dissolved  solids  throughout  the  Basin.  Other 
parameters  of  water  quality  can  be  included  in  the  model  with  little 
effort.  The  model  incorporates  hydrologic  and  water  quality  data  and 
the  Fiering-Pisaao  mathematical  model  described  in  the  report  "River 
Basin  Simulation  Program"  issued  by  the  Office  of  Comprehensive 
Planning  and  Programs,  March  1 9o 7 . 

Given  (1)  :he  River  Basin  Simulation  Program,  (2)  this  report, 
and  (3)  tape  of  operation- 1  hydrology,  other  investigators  can  study- 
various  combinations  of  v. ter  quality  management  schemes.  1 
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INTRODUCTION 


Purpose  and  Sccoa  of  Report 

This  report  describes  the  application  of  the  Fier ing-Pisano 
mathematical  model  to  the  Red  River  of  the  North  Basin.  The  report 
describes  pertinent  hydrologic  factors,  waste  loadings  to  the  river 
system,  and  methods  of  applying  the  model  to  determine  effects  of 
water  management  practices  in  the  Basin.  The  report  presents  the 
model  in  a  form  that  is  applicable  to  other  river  basin  systems  for 
which  appropriate  data  are  available. 

Need  to  Model  the  Red  River  of  the  North 

The  Red  River  of  the  North  and  its  tributaries  are  the  major  water 
supplies  for  over  160,000  people  and  several  industries  in  North  Dakota 
and  Minnesota,  In  return,  treated  and  untreated  v'astes  from  muni.cipalitic 
and  industries  ere  added  to  these  streams  at  various  points  in  the  Basin. 
Many  of  the  streams  frequently  contain  excessive  concentrations  of  dis¬ 
solved  salts  from  natural  sources.  Even  higher  concentrations  of  salts 
will  be  added  to  Sheycnne  River  and  Red  River  as  a  result  of  inflow  of 
planned  irrigation  return  flows  and  flushed  saline  waters  from  the 
Devil's  Lake  Basin. 

As  a  result  of  wastes  currently  added  Xo  the  system,  the  Red  River 
of  the  North  is  polluted.  Because  of  the  pollution  and  because  the  Red 
River  of  the  North  is  an  interstate  and  international  stream,  a  pollu¬ 
tion  abatement  enforcement  action  has  been  hold  and  the  river  is  under 
study  by  I '..VC  A  and  the  l.J.C.  (International  Joint  Commission).^ 


As  an  interstate  and  an  international  stream,  the  Red  River  of 
the  North  carries  the  wastes  of  dwellers  in  South  Dakota,  North  Dakota, 
Minnesota,  and  Canada.  At  the  same  time  it  represents  a  potential  of 
many  other  beneficial  uses  for  the  residents  in  the  Basin.  Water 
quality  standards,  currently  being  established  for  the  river,  must  be 
a  compromise  of  the  desirable  conditions  sought  by  the  users  of  the 
river  and  the  practical  limits  that  can  be  achieved  by  the  best  water 
management  practices.  Determination  of  the  values  representing  this 
compromise  can  be  accomplished  by  use  of  the  model  with  greater  speed 
end  validity  then  by  trial  end  error  or  by  other  time-consuming  detailed 
evaluations.  The  principal  purpose  of  the  model  is  to  explore  the 
impact  of  proposed  water  and  related  land  resource  developments  on 
water  quality  and  co  demonstrate  the  advantages  of  water  quality 
management  in  realizing  maximum  rccor.ivodations  of  conflicting  interests 
and  desires. 


WHAT  THE  MODEL  CAN  DO 


Attaining.  Water  Quality  Goals 

The  feasibility  of  attaining  some  set  of  instream  water  quality 
goals  can  be  investigated  in  the  following  fashion, 

Civen:  1.  a  set  of  water  quality  goals  that  may  vary 

spatially  and  seasonally  throughout  a  system. 

2.  a  configuration  of  economic  activity  which  is 
translatable  into  a  set  of  waste  loadings. 

3.  a  tool  that  (a)  preserves  the  essential  charac¬ 
teristics  of  the  river  basin  under  study  and 

(b)  is  capable  of  quickly  determining  the  conse¬ 
quences  of  imposing  the  above  loadings  on  the 
river  basin. 

Question:  Will  the  intended  goals  be  achieved  wi th  a  low 

probability  of  uneventful  failures?  If  not,  what 
blend  of  preventive  and  curative  pollution  control 
measures  must  be  initiated  such  that  the  frequency 
and  magnitude  of  violations  will  be  acceptable? 

Preventive  measures  may  require  either  a  curtailment  of  some 
desired  activity  or  the  installation  of  treatment  devices.  One 
curative  measure  that  could  be  applied  would  be  to  allow  pollution 
to  occur  and  then  to  correct  the  problem  by  either  scheduling  the 
transport  of  waste  such  that  maximum  advantage  is  taken  of  natural 
dilution  or,  alternatively,  by  scheduling  dilution  water. 
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Clearly  the  process  is  iterative.  A  change  in  goals  or  the 
introduction  of  a  new  blend  of  economic  activity  can  be  evaluated  to 
determine  whether  or  not  the  goals  will  be  met  without  pollution 
control  devices  or  to  determine  what  measures  must  be  taken  to  main¬ 
tain  the  goals.  The  different  sets  of  alternatives  nay  then  be  arrayed, 
condensed,  and  subsequently  reported  to  management. 

Management  then  can  determine  what  water  resources  management 
practices  should  be  impel emented  that  would  allow  North  Dakota  to 
develop  irrigation  and  flush  Devil’s  Lake  as  quickly  as  possible  while 
minimizing  the  number  and  degree  of  violations  of  any  proposed  water 
quality  TDS  (total  dissolved  solids)  standard. 

Sctiedul i  r.  g  of  Vlas  ta _D itch arge s 

The  consequences  of  imposing  any  seasonal  distribution  of  waste 
loadings  at  ten  different  sites  in  the  Basin  can  be  easily  evaluated. 
Although  the  model  is  limited  at  present  to  the  consideration  of  total 
dissolved  solids  and  instream  flow  requirements,  other  elements  can  be 
programmed  vi th  little  difficulty.  Maintenance  of  a  desired  level  of 
D.O.  (Dissolved  Oxygen) ,  for  example,  is  co.-m.only  expressed  as  a  certain 
flow  that  must  be  maintained.  Earlier  studies-^,  in  twin,  have  deter¬ 
mined  wk.'.t  minimum  flows  need  to  be  maintained  by  supplementing  natural 
flow  with  reservoir  releases.  Other  elements  can  be  reduced  to  cor., parable 
terms  and  evaluated. 

Flushing,  or  Devil’s  lake 

The  Bureau  of  Reclamation's  plan  is  to  use  irrigation  return  flows 
together  with  direcc  import  of  Garrison  Reset  toil'  water 

A 


to  raise  the 


levels  of  the  two  lakes,  thereby  diluting  the  salt  concentration  of 
the  lakes  to  an  acceptable  limit  for  the  placement  of  game  fish.  The 
time  horizon  necessary  to  accomplish  this  task  can  be  evaluated  by 
considering  Lho  effect  of  different  rates  of  import  from  the  two  major 
sources.  An  additional  consideration  may  be  that  of  storing  dilution 
waters  in  the  Lone tree  Reservoir  for  maintenance  of  water  quality  in 
the  Sheycnne.  Alternatively  or  jointly  different  operational  release 
schemes  for  the  existing  reservoirs  serving  the  Red  River  of  the  North 
may  be  tested  with  the  model. 

Management  of  Existing  Storage. 

Parameters  for  specifying  the  operating  rules  of  rule  curves  of 
existing  and  proposed  reservoirs  are  dr-.ta  input  to  the  model.  These 
rules  define  reservoir  operation  during  any  simulation  run.  Thus,  if 
irrigation  return  flo\:o  together  with  over flows  from  the  Devil's  Lake 
cause  serious  water  quality  problems,  or.c  course  of  the  investigation 
would  be  to  sec  if  existing  reservoirs  could  be  operated,  within  the 
bounds  of  their  existing  functions  and  in  a  slightly  different  fashion 
to  reduce  the  rather  of  violations  to  some  set  of  water  quality  goals. # 
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Tlic  Reel  River  of  the  North,  located  in  north  central  United  States 
and  sou tli  central  Cr.r.eGi,  begins  as  drainage  fro.a  South  Dakota ,  North 
Dakota,  and  Minnesota  at  the  junction  of  tier  Hois  dc  Sic-:::  and  Oliertail 
Rivers  at  Wnhpcton,  N.  Dak.,  and  Rrockcnr  id'  -  ,  Minn.  It  courses  400 
miles  northward  belvc-en  North  Dakota  and  Mf.nno.i-ota  throng’s  the  broad 
flat  plain  bed  of  the  glacial  lake,  Lake  Agassiz,  past  the  cities  of 
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record  (1912-1960)  in  1950  for  c:«.r..ple,  at  95,500  cfs  near  the  boundary, 

2 

represents  only  about  2.4  cfsm  runoff  --  considerably  lover 

2 

than  the  10  or  more  cfsm  runoff  common  to  the  more  hui.'.id  regions  of 
the  eastern  United  States. 

Highest  flows  coim.ionly  occur  during  the  spring,  in  March  and  April, 
from  snov:  and  ice  melt.  High  f  1  ovs  occur  also  during  the  nonvinter 
months  as  a  result  of  storms.  Flooding,  common  in  the  lower  reaches  as 
a  result  of  these  high  flows,  is  compounded  during  the  spring  periods  by 
the  progressively  northward  thawing  of  snow  and  ice  as  the  winter  ends. 

The  lower ,  northern  reach  of  the  river,  blocked  wi th  winter  ice  is 
incapable  of  conveying,  without  flooding,  the  large  flows  resulting 
from  melts  of  the  earlier  thawing  southern  sections  of  the  Basin. 

Low  flows  occur  frequently  i n  the  Basin  each  year,  generally  during 
late  fall  and  during  the  winter.  Seme,  st  all  streams,  unless  augmented 
by  reservoir  releases,  cease  tc>  flow  during  the  winter  months.  A  descrip¬ 
tion  of  flow  characteristics,  including  reletii.nf.hip*.  between  flow  and  IDS 
concentrations,  is  provided  for  several  streams  in  the.  section  describing 
the  mathc--.it  ical  model  of  the  Basin. 

Thirty  percent  of  the  average  flow  in  the  led  River  of  the  North  at 
Er arson,  Canada,  passes  through  control  structures  of  reservoirs  on  major 
tributaries.  Five  reservoirs  on  tributaries  in  the  Basin,  villi  a  total 
capacity  of  over  ?  million  acre-feet,  provide  some-  flood  relief.  Con¬ 
sideration  of  voter  quality  in  the  operation  of  these  reservoirs  could 
make  them  effective-  tools,  for  management  of  water  quality  in  the  riw  r . 
Table  I  provider,  information  on  these  reservoirs. 
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Storage  and  Operation  of  Five  Reservoirs 
in  the  Red  River  of  the  North  Basin 


RESERVOIR 

USABLE  STORAGE 

_ LLD. _ 

STORAGE  TI 
(YRS) 

ORWELL 

OTXEKTAIL  R. 

20, A 00 

0.1 

LAKE  TRAVERSE 

BOIS  DE  SIOUX  R. 

137,000 

2.5 

LAl'.E  ASHTABULA 
SHEYENKE  R. 

69,100 

0.9 

RED  LAKES 

RED  LAKE  R. 

1,905,000 

6.6 

hors 

SO.  BR.  PARK  R. 

3,550 

0.2 

OPERATING  SWISS 


ORWELL 

TRAVERSE 

ASHTABULA  Rl 

ED  LAKES 

KG'. 

. (cfO. 

OCT 

0.69) 

0.412 

0.595 

0.54  9 

300 

NOV 

.635 

.412 

.  580 

.545 

300 

DEC 

.550 

.412 

.  560 

.540 

300 

JAN 

.A  50 

.412 

.  540 

.  534 

300 

FEB 

.300 

.412 

.5)0 

.528 

300 

KVt’v 

.200 

.4  36 

.500 

.522 

300 

APR 

.156 

.900 

.595 

.545 

ISO 

RAY 

.  2.03 

.500 

.595 

.510 

ISO 

JUN 

.230 

.440 

.600 

.567 

3  80 

JUi. 

.  280 

.4  20 

.600 

.  564 

ISO 

ARC 

.420 

.41? 

.600 

.557 

1 80 

SEE 

.530 

.412 

.600 

.554 

ISO 

NOTE; 

Eigi 

LP.es  crc 

proportion  or 

tota.l  c.prnity  to  ! 

lined,  e..i' 

the: 

«  for  llo 

p'.r  e  Rer-orvoir , 

vhioh  arc  I'onthly 

drafts.. 

in  citVc 

per  *  c cord  . 
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Although  many  tributaries  contribute  to  the  flow  of  the  Red  River 
of  the  North,  those  from  Minnesota  have  a  greater  impact  on  the  flow 
of  the  river  than  those  from  North  Dakota.  Information  on  major  tribu¬ 
taries  is  provided  in  the  section  describing  the  model  of  the  Basin. 
Other  surface  water  features  of  the  Basin  include  many  small  lakes  in 
the  southeastern  section  of  the  Basin  that  drain  to  tributaries  of  the 
Red  River,  several  small  lakes  and  "pot  holes"  that  offer  no  direct 
surface  water  contribution  to  the  river  system,  and  the  closed  Devil's 
Lake  Basin  that,  at  present  is  not  continuous  with  the  river  system. 
Planned  water  resource  and  land  development  in  which  irrigation  return 
flows  will  be  routed  through  the  Devil's  Lake  Basin  into  the  Sheyenne 
River  will  have  a  major  impact  on  water  quality  of  the  Red  River. 

The  planned  process  of  routing  water  through  the  Devil's  Lake 
Basin  will  be  accomplished  by  diverting  about  750,000  acre-feet  of 
water  per  year  from  the  Missouri  River  to  irrigable  lands  north  of  the 
Basin.  Return  flows  from  this  irrigation,  although  expectedly  high  in 
mineral  content,  will  be  considerably  lower  in  concentrations  of  dis¬ 
solved  solids  than  water  in  the  Devil's  Lake  Basin  (concentrations  of 
TDS  generally  exceed  30,000  ppm),  and  are  expected  to  "freshen"  the 
water  in  the  Basin  lakes.  By  controlled  releases  of  the  freshened 
lake  waters  and  Garrison  Reservoir  water  into  Sheyenne  River  while 
continually  adding  water  to  the  lake  system,  the  lakes  arc  expected 
to  be  "freshened"  to  the  point  of  providing  a  beneficial  habitat  for 
fish  and  wildlife  without  seriously  impairing  water  quality  in  the 
Sheyenne,  and  eventually  the  Red  River  of  the  North. 
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Water  quality  problems  arc  not  limited  to  the  Devil's  Lake  area 
in  the  lied  River  Basin.  Concentrations  of  YDS  and  specific  ions  such 
as  chloride  and  sulfate,  frequently  exceed  desirable  limits  in  several 
streams  in  the  Basin.  lli-gb  concent  re' ion:;  of  salts  in  the  streams  arc 
attributable  primarily  to  the  mm  •  i .  1  composition  of  the  glacial  drifts 
and  sandstone  through  vhich  the  stre  a,  pass  or  from  which  they  receive 
groundwater  flow.  This  is  compounded  by  high  evaporation  rates  and  the; 
addition  of  municipal  and  industrial  wastes  to  the  streams. 

Water  Use 

Approximately  ?3  ir&d  (million  gallons  per  day)  of  water  was  used 

2/ 

for  municipal  supplies  in  the  Red  River  of  the  Worth  Basin  in  19 CA  S' 

Of  this,  about.  15  w.gd  we:,,  obtained  from  surface  water  sources;  the 
remainder  was  obtained  fror.  veils.  Water  for  industrial  use  amounts  to 
about  18  mgd  and  is  obtained  fro:-  veils  and  municipal  supplies.  An 
additional  amount  of  water,  about  AO  r-gd,  is  withdrawn  from  stre;  .s 
and  used  for  cooling  purposes  by  several  tli:  vro~clcc.tr  ic  plants. 


Poll  ell  on 

Red  River  of  the  north  and  way  of  its  tributaries  are  polluted. 
Wastes  fro. i  four  large  sugar  barf  wills,  several  potato  pi  o.  easing. 
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and  pathogenic  orga.ui  sm:: .  As  a  melt,  the  vatei  must  receive  a.ho-.v 
noi .  i.'.l  treat.  ent  before  it  is  used  fc”  :  .inicipnl  supplies  to  a::u’i 
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adequate  safeguards  to  health.  Also  because  of  the  pollution,  some 
stream  reaches  are  devoid  of  fish  and  the  v;atcr  constitutes  a  health 
hazard  to  persons  participating  in  water  contact  sports. 

One  complication  of  waste  disposal  in  the  Red  River  of  the  Korth 
system  is  the  long  winter  period  in  which  ice  covers  the  streams. 

Wastes  generally  are  held  in  lagoons  during  this  period  and  released 
during  the  high  runoff  period  of  the  spring  thaw.  Although  flows  are 
high  at  this  tire,  the  water  often  contains  high  concentrations  of 
dissolved  solids.  The  snow,  in  melting,  picks  up  and  transports 
saline  residues  from  the  land  surface  and  stream  banks  to  the  receiv¬ 
ing  streams. 

At  the  re rues t  of  the  Governor  of  Minnesota,  FWPCA  began  a  pollu¬ 
tion  study  of  the  Red  River  of  the  North  in  1964.  Efforts  are  being 
made  to  improve  water  quality  conditions  in  accordance  with  recoi.  •onda- 
tienn  of  the  conferences.  (Reference  to  Con for  dices ,  In  the  Matter  of 
Pollution  of  Interstate  Waters  of  the  Red  River  of  the  North,  Sept .  14 , 
1965.) 
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THE  MODEL 


Background 

A  very  general  anc’  flexible  package  entitled,  "River  basin 
Simulation  Program"  ~~  was  used  in  this  study.  It  is  a  scries  of 
programs  which  accept  certain  kinds  of  data.  The  package  end  the 
data  attempt  to  capture  the  pertinent  underlying  behavior  or  charac¬ 
teristics  of  the  process(os)  being  investigated.  Tin's  package,  and 
data  for  another  river  basin  would  again  become  the  model  of  that 
basin. 

Briefly,  the  inf  on  >al:ionnl  requi  rerents  for  the  Red  River  of  the 
North  model  include  data  on  (1)  historical  flown  at  gaging  stations; 

(2)  geometric  location  of  reservoirs,  waste  water  input  s,  and  voter 
users;  (3)  background  water  quality  relationships;  ov: novation  at: 
dansi.lcs;  (5)  r'.gnit. udes  of  water  u  c,  wane  input,  and  j  e.' .  r-.-oir 
volumes,  and  (6)  waste  scheduling  and  reservoir  rana.gc:  pract  ice  ••• . 

Ob  j  ec  t iy e 

The  objective  of  t  his  report  is  to  develop  a.  ra.U1  vlic’i  pi  cs.-t  vos 
the  stat  istical  characteristics  of  i.onthly  sire.  .  Per:  ar  ’  r.v.thly 
concentrations  cf  total  dissolved  solids  th>  ingli'Mi  th  l'a  •'  t  St 
portion  of  the  ltd  River  of  the  Not  I  h  syr  ten  c.o.  ’pa  i  t...  d  of  3,030 
of  river  divining  -'i0,000  squire  nilis.  The  nodel  i.ur.t  ft  .  .-id.-r  (!) 
natural  strea.  •flaws,  (2)  existing  reservoirs  cut  their  cement  ge¬ 

nu’.  (3)  existing  t 'r  .a1  s  a.n.l  waste  la  tern  flow  . 


1  2 
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meat  p: act  ices, 


Strat  cpy 


A.  Hydrologic 

1.  General 

The  hydrologic  model  was  defined  by  the  monthly 
flov.’  data  at  twelve  gages  located  throughout  the  Basin. 

Ten  of  the  gages  are  located  at  points  where  there  has 
been  no  flov;  regulation;  the  other  two  gages,  at  Ottertail 
River  below  Orv/ell  Dam,  near  Fergus  Falls,  Minnesota,  and 
at  Park  River  near  Grafton,  North  Dakota,  are  affected  by 
regulation.  Backrouting  was  performed  on  the  data  from 
the  latter  two  gages  to  restore  them  to  the  natural  regime. 

It  is  our  assumption  that  the  system  is  sufficiently 
described  by  (1)  using  the  twelve  gages  to  define  the 
Basin-wide  hydrology,  and  (?.)  superimposing  on  this  all 
flow  regulations  consisting  of  (a)  existing  reservoirs, 

(b)  existing  municipal  and  industrial  diversions  and 
return  flows. 

2.  Gages 

Shown  on  Figure  1  are  twelve  gages  considered  to  be 
(1)  sufficiently  representative  of  the.  various  distribu¬ 
tions  of  runoff  characteristic  in  the  Basin;  and  (2) 
adequately  widespread  as  to  capture  the  attendant 
regional  spatial  correlations  that  exist.  Shown  also 
on  Figure  1  are  tlr  areas  of  ungaged  portions  of  the 
Basin  that  were  represented  by  areal  transforms  (ndjuslcd 
in  some  instances  by  known  yield  i nf orrw: l i on) .  Roughly 
65  percent  of  the  total  area  was  not  gaged. 
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Record  length  between  years  1941-1960  was  chosen  to 

(1)  reduce  the  effort  and  the  associated  errors,  inherent 
in  record  extension  and  backrouting  of  regulated  flows; 

(2)  be  representative  of  vet  and  arid  periods,  and  (3) 
provide  a  period  when  each  of  the  twelve  gages,  save  one, 
recorded  information.  The  data  gap,  located  at  Coopcrstown, 
North  Dakota  was  filled  in  with  graphical  analysis  of  an 
adjacent  point.  The  mean  monthly  discharges  and  standard 
deviations  for  each  gage  are.  shown  on  Figure  2. 

Using  statistical  parameters  derived  from  the  monthly 
historical  flows  at  the  twelve  gages,  a  tape  of  two  thousand 
years  of  operational  hydrology  was  prepared  for  future  use. 
Pertinent  input,  and  output,  data  arc-  available  on  request. 

A  brief  description  of  the  technique  is  provided  in  the 
Appendix. 

3.  System 

Shown  in  Figure  3  is  the  system  of  components  used  in 
the  analysis.  Note  also  the  inset  table  which  describes: 

(1)  each  coordinate,  (2)  the.  supporting  drainage  area,  (3) 
the  gage  used  to  define  the  inflow  into  this  coordinate, 
and  (4)  the  factor  used  to  scale  the  flow  at  the.  gage.  It 
should  bo  clear  that,  the  50  coordinates  shown  were  more  than 
enough  to  define  the  system.  It  would  be  possible,  at  a 
later  time,  to  expand  the  current  system  to  incorporate 
more  waste  sites,  test  points,  and  damsites.  At  present, 
tremendous  flexibility  is  inherent  in  the  system. 
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SYSTEM  JitPfJKSENTATION 
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FIGURE:  3 


4.  Reservoirs  and  Rule  Curves 

Table  1  gives  pertinent  information  for  the  five 
existing  reservoirs  in  the  system  including  the  five 
sets  c:'  reservoir  operating  rules  used  in  live  analysis. 
The  s  i-r.sot.al  ruler  either  uaintr.in  a  certain  pool  level 
(expressed  as  percentage  cf  capacity)  or  alternatively 
stipulate  a  certain  draft  to  be  ii.sde  during  that  time 


pericc.  This  information  was  furnish'd  in  V.  S.  Arny 
Corps  of  ingi:  :ers  report.-^ 

Verification  of  the  hydraulic  Routing 
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comparison  -  observed  and  computed  mean  monthly 

DISCHARGES  AND  STANDARD  DEVIATIONS  W.  Y.  ISA!  - 1030 


FIGURE  A 


Example  2  -  (1951-1960)  All  reservoirs  were  in  operation 

during  tliis  period.  Shown  in 
Figures  5  and  6  are  the  compari¬ 
sons  of  the  observed  monthly  flows 
at  the  Grand  Tories  and  Emerson  gages 
with  those  fron  the  model  at  the 
sane  locations  and  for  the  same 
time  periods. 

Figure  7  presents  the  observed  mean 
monthly  flows  (1951-1960),  1  one 
standard  deviation,  and  their  re¬ 
spective  medal  c ova  ter  parts  for  the 
gages  ft  Fr  arson,  Grand  For  Its,  fad 
Fargo.  The  drainage  area  gaged  by 
the  three  stations  is  respectively 
40,200,  30,100,  rata  6,800  severe 
miles.  Table  11  shows  a  comparison 
of  the  ob.crved  average  discharge 
for  the  period  or  rc-co/d  with  that 
produced  by  the  model  for  a  number 
of  strategic  points  in  tlv  cy.tr-.’.. 

6 .  Cone  1  v  r  i  on 

The  agreement  bf  tween  estimated  and  observed  fire  dg; 
implies  that  the  translation  and  touting  a:  r v  pt  ice.r-  arc 

val id . 
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COMPARISON  -  033ERVED  AND  COMPUTED  MEAN  MONTHLY 
iSCHARGES  VV.  Y.  i95!  -  196 
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FIGURE 


CC.VPAPrSCK  -  OBSERVED  AND  COMPUTED  MEAN  MONTHLY 
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FIG  UHL  7 


Quality  -  Total  Dissolved  Solids 

1.  Natural  -  Background 

Presented  in  Figure  8  are  the  locations  of  six  water 
quality  monitoring  stations  f row  which  data  were  obtained 
to  define  average  monthly  concentration  --  discharge 
relationships.  Shown  also  are  the  areas  which  are 
assumed  to  exhibit  similar  leaching  and  discharge  behavior. 
The  pertinent  coefficients  for  these  relationships  are 
shown  in  Table  Ill. 

2.  Existing  waste  loadings 

Present  municipal  and  industrial  (potato  and  sugar  beet ) 
diversions.,  waste  return  flows  and  1  or.d i r.gs  for  the  bangs.::.: 
Fargo,  and  Grand  Forks  complexes  arc-  shown  in  Table  IV, 

The  pollution'll  cent  ' :  but  ?  c- •••  front  the  now  sugar  he  at  proces¬ 
sing  plant  at  Drryl.v.  was  not  considered. 

3.  Verifying  Quality 

For  the  years  H’:iC-lSoO,  eg.' a  3  i  tv  (as  TbS)  wcas-iae  cats 
ware  recorded  at  Lisbon,  Fo.-go,  a:  s'  Gruti;*  Folks,  New  th  ibdn.t 
The  i.  odd  was  operated  f  o:  the  s.  a  period  of  record.  She'...-, 
in  Figure  9  are  the  observer;  time  •••  weighted  r-o..thIy 
+  one  standard  deviation,  fi  the  system.  counterparts  fc>" 
Lisbon,  Fargo,  red  C.\  ad  Forks  .  In.  Figure  19  the  arc.--,  rl  • 
tive  probability  cm: res  f  uv  the  three  1 oea  t i ons  ere  gi\.n. 
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TABLE  II 

Comparison  of  rhc  Observed  and  Computed  Mean  Discharge 
for  Period,  Water  Years  1941-1960,  at  Selected  Gages 

•  MEAN  DISCHARGE 

STREAM  AND  LOCATION  Observed  Computed  Difference 


SsJ;il _ 

_ (cfs) 

(percent) 

Red 

River  of 

Che  North  at 

Wahpcton,  N.  Dak. 

480 

471 

-1.9 

Red 

River  of 

the  North  at 

Fargo,  N.  Dak. 

590 

578 

-2.0 

Red 

Lake  River  at  Crooks Ion,  Minn. 

937 

928 

-1.0 

Red 

River  of 

the  North  at 

Grand  Forks, N. Dak, 

2,193 

2,187 

-0.3 

Red 

River  of 

the  North  at 

Emerson,  Manitoba 

2,863 

2,887 

dO.S 

TABLE  III 

Coefficients  for  Relationships  of  Moan  Monthly  Concentration 
and  Discharge  for  Six  Monitoring  Stations 

(1> C--aQ^> T ,  vhere  C-concentration  of  total  dissolved  solids, 
in  ppm,  Q'di scharge ,  in  cfs.  L=-uppar  limit  of  concentration 

and  T-Iov:cr  limit  of  concentration.) 


Station 

a 

b 

L 

T  S: 

*.npl  c 

Breckeur i dge 

700 

-.1780 

600 

200 

50 

Cayuga 

1240 

-.2015 

2250 

300 

42 

Cooper st can 

985 

-.2155 

800 

200 

44 

Thief  R.  Falls 

950 

-.200 

1250 

2.50 

10 

Compos i to 

Snake,  Tamarac 
and  Tco  Rivers 

590 

-.170 

600 

160 

97 

Kalita  J  la 

890 

-.148 

760 

275 

77 
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LOG  AY  I  OH  MAP 

WATER  QUALITY  MONITORING  STATIONS 


1.  Vlilc1  Rico  River  n ;e.r  C?yc;ii 
c .  Ottvrleil  River  it  tree':;  r.ru'jc 

3.  Shoyc.iii.c  P.ivti  it  C c* c * j * c rstci.ai 

4.  Thief  River  el  Thief  River  (alls 

5.  Ccc’por.i  u>  Sr.ri.c,  levari-c ,  end  Tvo  Rivers 

6.  Pev!>i na  River  at  l.’i-ltnlle 


FIGU.M 


TABLE  IV 


Existing  Municipal  and  Industrial  Monthly  Diversions,  Waste  Return 
Flows,  and  Loadings  for  Wahpeton,  Fargo,  and  Grand  Forks,  N.Dak. 


WAHl'ETON  FARGO  GRAND  FORKS 

Diver-  Return  Concen-  Diver-  Return  Concen-  Diver-  Return  Conccn- 


Month 

si  on 
(cfs) 

Flo-..' 

(cfs) 

tration 

(ntft/1.) 

sion 

(cfs) .. 

Flow 

(cfs) 

tration 

(rue/1) 

sion 

(cfs) 

Flow 

(cfs) 

tration 
(mg/1 ) 

Oct . 

6.5 

A.  5 

800 

15 

9 

1000 

17 

15 

1000 

Nov. 

6.5 

4.5 

800 

15 

9 

1000 

17 

15 

1000 

Dec . 

3.0 

1.5 

800 

15 

4.5 

1000 

17 

15 

1000 

Jan . 

3.0 

1,5 

800 

15 

4.5 

1000 

17 

15 

1000 

Feb. 

3.0 

1.5 

800 

15 

4.5 

1000 

5 

5 

1000 

Mar . 

3.0 

1.5 

800 

9 

4.5 

1000 

5 

5 

1000 

Apr . 

3.0 

9.5 

800 

9 

30.0 

1000 

5 

5 

1000 

May 

3.0 

9.5 

800 

9 

30.0 

1000 

5 

5 

1000 

June 

3.0 

3.0 

800 

9 

9 

1000 

3 

3 

1000 

July 

3.0 

3.0 

800 

9 

9 

1000 

3 

3 

1000 

Aug. 

6.5 

4.5 

800 

9 

9 

1000 

3 

3 

1000 

Sept . 

6.5 

4.5 

800 

9 

9 

1000  ■ 

3 

3 

1000 

NOTH : 

An  ca 
Red  R 

t  i:.s  Tiled 
Ivor  of 

5?5  tons/year  of 
the  North,  sssi'n 

lir.e 
ed  cl  i  s 

(sugar  bee 
tribute J  a 

t)  is  di 
s  follow 

s charge 

d  into  t! 

Month 

Tons. 

Percent 

Month 

Tons 

Percent 

Month 

Tons 

Percent. 

Oc  t  . 

63.? 

13.0 

Feb. 

0 

0 

June 

105 

20 

Nov . 

18.4 

3.5 

Mur . 

105 

?  0 

July 

0 

0 

Dlc  . 

18.4 

3.5 

Apr , 

105 

20 

Aug . 

0 

0 

Jan . 

0 

0 

May 

105 

20 

Sept . 

0 

0 
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COMPARISON  -  OL' SERVED  AND  COMPUTED  MEAN 
MONTHLY  CONCENTRATION  Or  TOTAL  DISSOLVED 
SOLIDS  AND  STANDARD  DEVIATIONS  -V7.Y.  ID5S-TSSO 


f  lGUfCC  9 


COMPARISON  -  OBSERVED  AND  COMPUTED  FREQUENCY 
CURVES  OF  CONCENTRATION  OF  TOTAL  DISSOLVED  SOLIDS 
,  W.  Y.  I95G  -  I960 
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2  '0  20  iO  40  53  60  70  63  90  59 


FIGURE  10 
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appendix 


Mathematical  Model  of  Streamflow 
in  a  River  Basin 

One  scheme  for  a  mathematical  model  of  streamflow  in  a  river 
basin  is  used  in  a  FWPCA  computer  program,  "River  Basin  Simulation 
Program."  The  program  produces  synthetic  discharges  that  duplicate 
specific  characteristics  of  historical  streamflow  data.  A  long  synthetic 
record  can  provide  a  large  variety  of  possible  combinations,  sequences  of 
events,  extremes,  and  conditions  not  yet  documented.  The  long  record 
enables  the  user  to  define  a  wider  spectrum  of  discharge  events  than 
the  historical  data  provide,  to  test  a  design,  plan  a  facility,  schedule 
an  operational  pattern,  or  otherwise  consider  the  possible  range  of 
conditions  at  selected  locations. 

In  addition  to  the  sample  provided  by  the  historical  record,  addi¬ 
tional  samples  can  be  synthesized  thereby  permitting  the  user  to  better 
gage  the  reliabilities  of  various  alternative  courses  of  action.  In 
other  words,  one  can  (1)  measure  various  consequences  for  the  same 
course  of  action,  (2)  study  the  distribution  or  occurrence  frequencies 
of  the  consequences,  and  then  (3)  decide  on  the  desirability  of  the 
course  of  action. 

Three  specific  characteristics  of  historical  streamflow  data  are 
used  in  this  model  for  each  gaging  station  and  for  each  time  interval. 

The  time  interval  (season)  can  be  whatever  is  compatible  with  the 
particular  study  --  monthly,  bimonthly,  trimonthly,  semiannual,  or 
annual.  Calendar  months  are  usually  used  in  river  basin  studies. 
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Hie  first  parameter  is  the  ar  ithr.ietic  mean ,  in  any  units,  for 
each  season  for  the  period  of  historical  record  common  to  all  records 
in  the  basin.  Thus,  for  each  gaged  point,  a  set  of  twelve  mean  dis¬ 
charges  represents  the  central  value  of  the  distribution  for  each 
calendar  month. 

The  second  parameter  is  the  standard  deviation,  in  the  same  units 
as  the  mean,  for  each  season  which  is  taken  to  5*  the  square  root  of 
the  mean  of  the  squares  of  deviations  from  the  mean  of  all  discharges 
for  each  season.  This  parameter  provides  a  numerical  index  of  dispersion 
of  discharge  magnitudes  representing  each  of  the  twelve  months. 

The  third  parameter  is  the  lag  correlation  coefficient.  Because 
each  monthly  discharge  is  dependent  to  some  extent,  however  small,  on 
the  discharge  during  the  preceding  month,  the  discharges  are  not 
completely  random.  The  lag  correlation  coefficient  relates  the  serial 
interdependence  between  each  month  and  the  month  immediately  preceding. 
Thus,  a  set  of  12  lag  correlation  coefficients  for  each  gaging  station 
record  measures  the  temporal  (sequential)  relationships  of  the  historical 
data . 

The  hydrology  of  a  large  river  basin  varies  with  location.  Some 
influences  causing  the  variation  are  obvious,  others  are  subtle.  The 
transition  from  place  to  place  may  be  sudden  or  gradual .  Each  gaging 
station  within  a  whole  basin  is  a  sampling  point  ol  strer.mf low  compo¬ 
nents  of  the  hydrology  of  part  of  the  basin.  Gages  close  to  each  other 
tend  to  exhibit  similar  characteristics  and  arc  therefore  not  completely 
Independent  among  themselves.  To  consider  this  tendency,  the  model 
estimates  and  uses  a  table  (matrix)  of  correlations  of  each  gage  with 
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each  of  the  other  gages  for  each  of  the  twelve  months.  These  12 
tables  measure  the  spatial  relationships  among  the  gaged  discharges 
with  their  variation  from  month  to  month. 

The  temporal  and  spatial  parameters  derived  from  the  historical 
data  are  used  to  generate  a  multivariate  distribution  function  which 
in  turn  is  used  to  generate  synthetic  discharges  for  each  gaged  point 
in  the  basin.  Any  of  three  distributions  --  normal,  log  normal,  or 
gamma  --  can  be  selected.  The  computer  then  generates  a  large  number 
of  synthetic  discharges  (operational  hydrology)  based  on  the  assumed 
distribution  and  having  the  same  parametric  characteristics  (mean, 
standard  deviation,  lag  correlation,  and  spatial  correlation)  as  the 
historical  data.  The  parameters  of  the  synthetic  record  are  computed 
to  provide  a  comparison  with  those  of  the  historical  data  to  insure 
proper  replication. 


